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Abstract 

We give a short review of a number of different models for treating magnetization effects in plasmas. 
In particular, the transition between kinetic models and fluid models is discussed. We also give examples 
of applications of such theories. Some future aspects are discussed. 

1 Introduction 

The field of quantum plasmas has been rapidly growing over the last decade. In particular, studies regarding 
the nonlinear properties of systems in which quantum and collective effects play an important role have 
been in focus (for a review, see HI). However, there are also numerous studies of quantum plasmas where 
magnetization effects are of interest. Here, the intrinsic magnetic moment of the plasma constituents 
give rise to new collective dynamical properties due to indirect spin interactions (through effective field 
excitations) as well as spin-velocity couplings. 

The above physical systems can be described in a multitude of ways. Here we will give a short overview 
of part of such descriptions. We start with the "heuristic" approach of Madelung, for which a decomposi- 
tion of the system wave function into phase and amplitude leads to the definition of macroscopic density, 
velocity, and spin variables. We then go on to describe the more detailed effective field quantum kinetic 
theory, through which the relevant fluid moments may be defined and the concomitant fluid equations de- 
rived, as well as giving the opportunity to analyse proper kinetic effects in quantum plasma systems. We 
give a brief account of possible applications and results of the quantum fluid/kinetic models. 

2 The Madelung Approach to Quantum Dynamics 

A rather generic approach to quantum fluids is the use of a Madelung decomposition of the system wave 
function, in which the amplitude is translated into a density and the gradient of the phase determines the 
velocity variable. Such a decomposition will below be reviewed, and the results obtained will be compared 
later with the moment hierarchy obtained through a more rigorous quantum kinetics approach. 

2.1 The Schrodinger equation 

The basic equation of nonrelativistic quantum mechanics is the Schrodinger equation. The dynamics of an 
electron, represented by its wave function y/, in an external electromagnetic potential is governed by 

+ — vV+e^ = 0, (1) 
at 2m e 

where h is Planck's constant, m e is the electron mass, and e is the magnitude of the electron charge. This 
complex equation may be written as two real equations, writing y/ = y/n expiS/h, where n is the amplitude 
and S the phase of the wave function, respectively [2 |. Such a decomposition was presented by de Broglie 
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and Bohm in order to understand the dynamics of the electron wave packet in terms of classical variables. 
Using this decomposition in Eq. (Q]), we obtain 



dn 
dt 



-nV-v, (2) 



and 

m e — = eV$ + — V — ^ ), (3) 
dt 2m e \ y/n J 

where the velocity is defined by v = VS/m e and d /dt = d t + v • V. The last term of Eq. (0 is the gradient 
of the Bohm-de Broglie potential, and is due to the effect of wave function spreading, giving rise to 
a dispersive-like term. We also note the striking resemblance of Eqs. (f2]i and <[3j to the classical fluid 
equations. 



2.2 The Pauli equation 

The non-relativistic evolution of spin-1/2 particles, as described by the two-component spinor f( a ), is given 
by the Pauli equation (see, e.g., O) 



dt 



h 2 /„ ie ^ 2 



— ( V + -A 1 -/i B B-cr + c0 



Y = 0, (4) 



where A is the vector potential, jj.B = eh/2m e is the Bohr magneton, and a = (<7i , 02, 03) is the Pauli spin 
vector. 

Now, in the same way as in the Schrtidinger case, we may decompose the electron wave function y/ 
into its amplitude and phase. However, as the electron has spin, the wave function is now represented by 
a 2-spinor instead of a c-number. Thus, we may use y/ = y/ri exp(iS/ft)(p, where <p, normalized such that 
<p T <p = 1, now gives the spin part of the wave function. Multiplying the Pauli equation © by \j/ *, inserting 
the above wave function decomposition and taking the gradient of the resulting phase evolution equation, 
we obtain the conservation equations 

dn „ 

— = -nV ■ v (5) 



dt 



and 



dvi e . n 2 d {V 2 ^\ 

dt m e Lnr e axj \ y/n J 

HB s dBj h 2 d (6) 

m e 1 dxi 4m 2 n dxj 

respectively, where e,-,-^ is the fully antisymmetric (pseudo-)tensor and we have used Einstein's sum con- 
vention so that a sum over indices occurring twice in a term is implied, where i,j,k, ■■■ — 1,2,3. The spin 
contribution to Eq. (O is consistent with the results of Ref. Q3). Here the velocity is defined by 

1 eA 
v= — (VS-ih(p J V(p) + , (7) 

m e m e c 

the spin density vector is 

s = <p f (7<p, (8) 

which is normalized according tcQ 

|s| = l, (9) 
and we have defined the symmetric gradient spin tensor 

r, = ^. do) 

dxj dxj 



'An alternative choice for the normalization would be to choose |s| = ft/2. 
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Moreover, contracting Eq. <j4j by y^O, we obtain the spin evolution equation 



dsj ( 2/ig h 



dt \ h 7 2m e n 



d ( d 



dx[ y dxi 



s k - (11) 



We note that the last equation allows for the introduction of an effective magnetic field B e ff = (2/ig /fi)B — 
h/ [dj(ndjS)] I (2m e n). However, this will not pursued further here (for a discussion, see Ref. J2]). 

Comparing the effects due to spin from the Pauli dynamics with the Schrodinger theory, we see a 
significant increase in the complexity of the fluid like equations due the presence of spin. The fact that the 
spin couples linearly to the magnetic field makes the dynamical aspects of such Pauli systems very rich. 
Moreover, when going over to the collective regime, the back reaction through Maxwell's equation can 
yield interesting new properties of such spin plasmas. In fact, the introduction of an intrinsic magnetization 
can give rise to linear instability regimes, much like the Jeans instability. 

2.3 Collective plasma dynamics 

As pointed out in the previous section, the route from single wavefunction dynamics to collective effects 
introduces a new complexity into the system. At the classical level, the ordinary pressure is such an effect. 
In the quantum case, a similar term, based on the thermal distribution of spins, will be introduced. 

2.3.1 Multistream model 

The multistream model of classical plasmas was successfully introduced by Dawson 0. Here we will focus 
on the electrostatic interaction between a multistream quantum plasma described within the Schrtidiner 
model, a system first investigated in Ref. J5) (where also the stationary regime was probed). Thus, we have 
the governing equations (0 and (O but for N beams of electrons on a stationary ion background, i.e., Using 
this decomposition in Eq. (fTJ, we obtain 

— - = -« a V-v C( , (12) 



and 



dt 



d\ a e „ , h 2 



n 



a 



-JT = — Vtf> + ^V ), (13) 

dt m e lint 



now coupled through the self-consistent electrostatic potential governed by 

e N 

V 2 4> = - J>a-*>). (14) 

£ o a =i 

Here «o is the density of the stationary ion background. 

In the one-stream case (a = 1), we have the equilibrium solution v = vo (a constant drift relative 
the stationary ion background) and the constant electron density n = no (such that = 0). Perturbing 
this system a Fourier decomposing the perturbations, such that n = «o + 5«exp[;(k ■ x — cot)]), v = vo + 
5vexp[/(k- x — of)], and = 50exp[/(k - x— Of)], we obtain JUS] 

(o-k.v ) 2 = o^ + Sf, ( 15 ) 

where the last term is the Bohm-de Broglie correction to the dispersion relation. Here we have the electron 
plasma frequency o p = (e 2 /?o/£() OT e) 1//2 - 

Similarly to the one-stream case, we obtain the dispersion relation ||5]|6) 



(©-voi -k) 2 -^^4/4 



mt 



(o-vo2-k) 2 -ft 2 /fc 4 /4fi 



(16) 
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for two propagating electron beams (with velocities Voi and V02) with background densities «oi an d "02- The 
quantum effect has a subtle influence on the stability of the perturbed plasma. For the case «oi = «02 = "o/2 
and voi = — V02 = Vo, we have the instability condition 

w{ l -w) <H2< ^ (17) 

in terms of the normalized wavenumber K = kvo / (O p and the quantum parameter H = h(Q p /m e \>Q (see Fig. 
1) GJ 0. We see that when H = 0, we have unstable perturbations for < K < 1, but when H 7^ a 
considerably more complex instability region develops. 

A model for treating partial coherence in such systems, based on the Wigner transform technique QE 
l9l [l0l . can also be developed (51 (see also Ref. ifTTI ). Moreover, using the equations §5^ and ©, a similar 
framework may be set up for electron streams with spin properties. 



2.4 Fluid model 

2.4.1 Plasmas based on the Schrodinger model 

Suppose that we have N electron wavefunctions, and that the total system wave function can be described 
by the factorization y(xi i x 2, • • - x n) = Wi V2 ■ • • Wn- F° r each wave function i// a , we have a corresponding 
probability 2? a . From this, we first define y/ a = n a exp(iS a /Ti) and follow the steps leading to Eqs. (f2]i and 
(0. We now have such equations the wave functions { l/a}- Defining lfT2l 

N 

n=Y,^an a (18) 

a=\ 

and 

V^(Va)=L " > ( 19 > 

a=\ " 

we can define the deviation from the mean flow according to 

w«=v a -v. (20) 

Taking the average, as defined by $1% , of Eqs. (f2]) and (01 and using the above variables, we obtain the 
quantum fluid equation 

dn 
dt 

and 



-nV-v (21) 



d\ e „, 1 „ h 2 „ //V 2 A /^ a 



dt m e nm e 2m 2 , 



= - - — Vp + — - V , (22) 



a 



where we have assumed that the average produces an isotropic pressure p = m e n(\vr a \ 2 ) We note that the 
above equations still contain an explicit sum over the electron wave functions. For typical scale lengths 
larger than the Fermi wavelength kf, we may approximate the last term by the Bohm-de Broglie potential 

G2 



(23) 



Using a classical or quantum model for the pressure term, we finally have a quantum fluid system of 
equations. For a self-consistent potential we furthermore have 



V 2 = —{n-m). (24) 
£0 
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2.4.2 Spin plasmas 



The collective dynamics of electrons with spin and some of the spin modifications of the classical dispersion 
relation was presented in Ref. Ifl3l . Here we will follow Refs. Ifl3l and |[T4l for the derivation of the 
governing equations. Suppose that we have N wave functions for the electrons with magnetic moment 
fi e = —jJ.B, and that, as in the case of the Schrtidinger description, the total system wave function can be 
described by the factorization tyA = y/i y/2 ■ • • Wn- Then the density is defined as in Eq. ([T8l and the average 
fluid velocity defined by (T% . However, we now have one further fluid variable, the spin vector, and 
accordingly we let S = (s«). From this we can define the microscopic microscopic spin density y a = 
s a - S, such that (y a ) = 0. 

Taking the ensemble average of Eqs. © we obtain the continuity equation (fJTJ, while we the the 
ensemble average applied to © yield 



d\'i 
~dt 



{Ei + Eij k VjB k ) 



1 dp 

nm e dx; 



and the average of Eq. (fTTT i gives 



~dt 



— - —r-£ijkBjS k - 



| fi 2 d /v 2 yn\ 

1 dZij 



1 



-F 



spin 



spin 



(25) 



(26) 



respectively. Here the force density due to the electron spin is 



-.spin 



dBj Tr d r , . N1 



Am e dx: 



dSk\ /y ffi A | Jd&ak 
dxi J \ dxj I \ dxi 



dSk 
dxj 



consistent with the results in Ref. IfTSI , while the asymmetric thermal-spin coupling is 

£y = nm e (.y ai w aj ) 

and the nonlinear spin fluid correction is 



(27) 



(28) 



n spm 



; £ ijkSj 



nh 



-£ijk 



t dSk 
dx a \ dx a 
d 
dx a 



-£ijkSj 



9 ( n /dy ak 

dx u \ \ dx a 



aj 



n a ^—{Sk + y a k) 
dx a 



(29) 



where Tij = (djS a )(djS a ) is the nonlinear spin correction to the classical momentum equation, fy = 
( (^i^(a) a ) (djy^ )) is a pressure like spin term (which may be decomposed into trace-free part and trace). 
We note that, apart from the additional spin density evolution equation 1261 . the momentum conservation 
equation dZSb is considerably more complicated compared to the Schrtidinger case represented by ( [22} . 
Moreover, Eqs. d25b and ( 126b still contains the explicit sum over the states, and has to be approximated 
using insights from quantum kinetic theory or some effective theory. 

The coupling between the quantum plasma species is mediated by the electromagnetic field. By defini- 
tion, we let H = B//io — M where M = —ln^gS/h is the magnetization due to the spin sources. Ampere's 
law V x H = j + £q^E takes the form 



1 <9E 

VxB = ^ (j + VxM) + - — , 
c L at 

where j is the free current contribution The system is closed by Faraday's law 

VxE = —. 

dt 



(30) 



(31) 
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2.5 The magnetohydrodynamic limit 

The concept of a magnetoplasma was first introduced in the pioneering work [161 by Alfven, who showed 
the existence of waves in magnetized plasmas. Since then, magnetohydrodynamics (MHD) has found 
applications in a vast range of fields, from solar physics and astrophysical dynamos, to fusion plasmas and 
dusty laboratory plasmas. 

Magnetic fields, an essential component in the MHD description of plasmas, also couples directly to 
the spin of the electron. Thus, the presence of spin alters the single electron dynamics, introducing a 
correction to the Lorentz force term. Indeed, from the experimental perspective, a certain interest has 
been directed towards the relation of spin properties to the classical theory of motion (see, e.g., Refs. 
H7] H8] HU EH! EH E2] El] EU ESI EU Ell EU |29]). In particular, the effects of strong fields on single 
particles with spin has attracted experimental interest in the laser community |fl9l l20l |2"T1 l22l l23l l25l . 
However, the main objective of these studies was single particle dynamics, relevant for dilute laboratory 
systems, whereas our focus will be on collective effects. 

We will now include if the ion species, which are assumed to be described by the classical equations 
and have charge Ze, we may derive a set of one- fluid equations |[T4ll . The ion equations read 

dm „ 

-i = - ni W-M, (32) 
dt 

and 

mm ^ = Zem (E t + e ijk v Ij B k ) - , (33) 

Next we define the total mass density p = (m e n + mjnj), the centre-of-mass fluid flow velocity V = 
(m e nv e + mjnj\j)/ p, and the current density j = — en\ e + Zeri]\j. Using these denfinitions, we imme- 
diately obtain 

f = -pV-V, (34) 

from Eqs. (f2~TT > and (132) . Assuming quasi-neutrality, i.e. n w Zn/, the momentum conservation equations 
(l25l > and (J33j give 

dV, dUij dp Zh 2 P d /V 2 VP\ spin 

P^7 - £ ijkJj B k~- -3 3— + 7. 3— — , (35) 

at axj axi Zm e mj oxj \ ^Jp ) 

where 11 is the tracefree pressure tensor in the centre-of-mass frame, and P is the scalar pressure in the 
centre-of-mass frame. We also note that due to quasi-neutrality, we have n e w Zp /mi and v = V — m/j/Zep, 
and we can thus express the quantum terms in terms of the total mass density p, the centre-of-mass fluid 
velocity V, and the current j. With this, the spin transport equation (TZoT i reads 

dS, = m.j d* ?^ _ m d^l ^ nspin 

H dt Ze Jj dxj ti Z dxj Z 

In the momentum equation ( 1351 ), neglecting the pressure and the Bohm-de Broglie potential for the sake 
of clarity, we have the force density j x B + F sp ; n . In general, for a magnetized medium with magnetization 
density M, Ampere's law gives the free current in a finite volume V according to 

j = — VxB-VxM, (37) 
Mo 

where we have neglected the displacement current. The surface current is an important part of the total 
current when we are interested in the forces on a finite volume, as was demonstrated in Ref. |[T4ll and will 
be shown below. 

It it worth noting that the expression of the force density in the momentum conservation equation can, 
to lowest order in the spin, be derived on general macroscopic grounds. Formally, the total force density 
on a volume element V is defined as F = limv^o{Y.afoc/V), where f a are the different forces acting on 



6 



the volume element, and might include surface forces as well. For magnetized matter, the total force on an 
element of volume V is then 



v 



{Mxn)> 



where (neglecting the displacement current) j to t = j+ VxM. Inserting the expression for the total current 
into the volume integral and using the divergence theorem on the surface integral, we obtain the force 
density 

F tot =jxB+M*VB*, (39) 

identical to the lowest order description from the Pauli equation (see Eq. (l35ll). Inserting the free current 
expression d37l i. due to Ampere's law, we can write the total force density according to 

B \ — '-rink-. 



F' = -d'[ M • B + d k (H'B k ). (40) 

V 2 Mo / 

The first gradient term in Eq. d40l > can be interpreted as the force due to a potential (the energy of the mag- 
netic field and the magnetization vector in that field), while the second divergence term is the anisotropic 
magnetic pressure effect. Noting that the spatial part of the stress tensor takes the form [TT31 

T ik = —H i B k + (B 2 /2hq - M ■ B)S ik , (41) 

we see that the total force density on the magnetized fluid element can be written F' = —d/ ( T' k , as ex- 
pected. Thus, the Pauli theory results in the same type of conservation laws as the macroscopic theory. The 
momentum conservation equation ( f33T > then reads 

p ( ^- + V • v) V = -V (^— -M-b] + B k d k H - Vp, (42) 



dt J \2no 

where for the sake of clarity we have assumed an isotropic pressure, dropped the displacement current term 
in accordance with the nonrelativistic assumption, and neglected the Bohm potential (these terms can of 
course simply be added to d42b). This concludes the discussion of the spin-MHD plasma case. Next, we 
will look at some applications of the derived equations. However, it should be noted that in many cases the 
spins are close to thermodynamic equilibrium, and we can thus write the paramagnetic electron response 
in terms of the magnetization lTT4l 

M^tanhf^B, (43) 
mi \k B T J 

instead of using the full spin dynamics. Here B denotes the magnitude of the magnetic field and B is a unit 
vector in the direction of the magnetic field, Icb is Boltzmann's constant, and T is the electron temperature. 



3 Spin Quantum Kinetics 

Quantum mechanics in terms of quasi-distribution functions is perhaps the formulation with the closest 
resemblance to classical statistical mechanics. The formulation started with Wigner's paper [7 | together 
with Weyls correspondence principle 1301 . In terms of this formulation the state of the system is no longer 
described by a density matrix but instead a phase-space distribution function. Similarly operators are trans- 
lated into phase functions. Calculating the expectation value of an operator is then a matter of calculating a 
phase space integral a corresponding function weighted by the distribution function. The method has been 
applied to a wide range of problems. For example it has been applied in optics lF3Tl[32ll . collision theory 
1 33 34], nonlinear theory |35l and transport problems in solid state physics, see for example |36l [371 and 
references therein. 

There are many different ways to define a quasi distribution function in quantum mechanics. The most 
known examples are probably the Glauber-Sudarshan p-distribution [38 39], the q-distribution (40j and 
the related Husimi distribution RD . The many different definitions basically comes from the fact that 
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the position and momentum operators do not commute, so the transformation between an operator and a 
phase space function is not unique. See for example l42l for a review of the different phase space distri- 
bution functions. The Wigner distribution corresponds to ordering the position and canonical momentum 
operators symmetrically. 

When considering a particle in a magnetic field gauge invariance has to be assured. This can be done by 
adding a phase factor to the definition of the Wigner function J43). The phase factor will then compensate 
for the change of phase that occurs in the density matrix when performing a gauge transformation. When 
dealing with the gauge invariant Wigner-Stratonovich the Weyl correspondence is modified l44l . The 
natural variables to use are the position and kinetic momentum. 

The formulation of quantum mechanics in terms of phase space distributions has also been generalized 
for spin particles |45 [ 451 |47l [48l . Also in this case the definition is not unique and one can find analogs to 
the different definitions in the phase space case B31 . 



3.1 Scalar quasi-distribution theory for a spin plasma 

The distribution function which we will work with here is the combination of the Wigner distribution 
function [7J for the phase-space variables and the q-function for the spin degree of freedom l47l . This 
combination of distribution functions was used in |49l and it turns out to yield a intuitive description of 
spin- 1/2 particles in an extended phase space. 

Given a 2-by-2 density matrix (in the Schrodinger pricture p(x,y,f) for a spin particle the extended 
phase-space distribution function is defined by 

/(x,p,s,0 = ^Tr[(l + s- <T)W(x,p,0] , (44) 

where Tr denotes that the trace is to be calculated of the resulting 2-by-2 matrix and where cr is a vector 
with the Pauli matrices as components, s is a vector on the unit sphere. The Wigner distribution matrix 
function is given by 

where q is the charge of the particle and A is the vector potential. The integral over the vector potential 
is there to ensure gauge invariance. The momentum variable p is the gauge invariant kinetic momentum 
related to the canonical momentum p c by p = p ( — qA(x,t). This distribution function is defined on an 
extended phase-space (x, p, s) and can in principle be used to calculate the expectation value of any observ- 
able defined by an operator O = 0{k,p, cr). The way to do this is to use the modified Weyl correspondence 
B4l together with a transformation for the spin variable to obtain a phase-space function, and subsequently 
take the average of this function weighted by the distribution function. For an operator depending on x,p 
the corresponding phase-space function is obtained by 



Ofep) = ffyeMntf***™*] ( x + 1| 6 |x- |) 



(46) 



It can sometimes be found more easily by first putting the position and kinetic momentum operators of the 
operator O in symmetric order using the commutation relation and then make the substitution x — >• x and 
p — > p. For example the pressure tensor which just contains the kinetic momentum operator is obtained 
by Pjj = pipj — s- piPj. Note that the gauge dependent Wigner distribution this correspondence is not so 
simple anymore, since we are then dealing with the position and canonical momentum operators which has 
to be ordered symmetrically using the commutation relation. So for the moment above, for example we 
have to put the combination Py = \pd — qAj(x)][p C j — qAj(x)] in symmetric ordering and then make the 
substitution i, — > a, and /3,- — > p,-. This is a quite difficult task for a general operator since the form of the 
vector potential is not necessarily known. 

The spin space function corresponding an operator O depending on cr is obtained by transformation 

°( s )= 4 V rr [ (i+s '' T) ^- (47) 
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If the operator in question depends on both the position and momentum operators and the spin, both of 
these transformations have to be made, see Ref. []. Since er 2 = 1 the only possible spin operators we may 
have is the identity operator 1 and er. The transformation (|47T > above yields 

6 = er -> 0(s) = 3s. (48) 

The momentum variable in the Wigner function above is the canonical momentum. In the presence of 
a magnetic field it is often more convenient to work with the gauge invariant kinetic momentum. 

3.2 Evolution and the long scale length limit 

The Hamiltonian for a spin- 1/2 particle in a magnetic field is given by 

Lm 

where V and A are the electromagnetic potentials, B = V x A is the magnetic induction and fj,g is the 
magnetic moment of the particle. Specifically, for an electron, the magnetic moment is given by /J-b = 
qgh/(4m) where g ss 2.001 is a correction factor deduced from quantum electrodynamics l50l . Note that 
we have used q = — \e\ so that the magnetic moment is negative. The evolution equation for the density 
operator is given by the von Neumann equation 

ihd,p = [p,H}. (50) 

Taking the transform d44b of this equation, it is possible to derive an evolution equation for the extended 
phase space Wigner function /(x,p,s,f ), see Ref. fl49l and it is given by 

i + (v+Af)-V^+i[E+(v + v)xB)].V v / 
dt m 

+— V,[(s + V,).fi].V v /+^(sxfi).V,/ = 0. (51) 
m n 

where v = p/m is the velocity and we have defined the operators 

E = E(x) / rfscos — V. v -V v (52) 
J -1/2 \m J 

B = B x) / djcos — V x -Vv (53) 
J -1/2 \m ' J 

^-B(x) [ l/2 dsssmf—^Vx-^^) -V v . 



Av = -S-B(x) / dss sin — V r -V v -V v . (54) 

TH J -1/2 \ OT / 

Note the similarity of the equation above with the classical Vlasov equation. In order to compare it further 
we may consider the semiclassical limit which is applicable when the typical length scale is much shorter 
than the de Broglie wave length Ti/ (mv), where v is the typical velocity of the system. We may then expand 
the sine and cosine operators above. Keeping terms up to order fi the resulting equation is 

TP + v • V x f + — [E + v x B] + — V r [(s + V,) • B] • V v / + ^ (s x B) • V. v / = 0. (55) 
dt m m n 

The first term proportional to the spin s is the dipole force of a magnetic moment in an inhomogeneous 
magnetic field and the last term is accounts for the spin precession. Both of these can be understood from a 
classical analog, however, there is also a quantum dipole term given by (fls/m)'V x (B ■ V s .) • V,,/. This term 
accounts for the fact that the spin is not a classical dipole and we can for example not have a distribution 
function proportional to 8 (s — z), where £ is a unit vector in the z-direction. The corresponding distribution 
function in the quantum mechanical case is proportional to 1 + cos 9 S where 9 S is the angle between the 
spin direction and the z-axis. The macroscopic magnetization the classical and quantum cases are still the 
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same which is ensured by the factor 3 in the transformation ((48) above. This factor occurs in the quantum 



M q (x)=3MB / £/£2/(x,v,s) 



but is absent in the corresponding classical equation 



Md(x) =Hb / dfls/ d (x,v,s), 



(56) 



(57) 



where f c \ is the classical distribution function and we in both cases have d£l = d 3 vd 2 s. 

A semi classical version of Eq. ( f55l l (where the quantum dipole term is missing) has been applied to 
find a new type of resonance due to radiative corrections to the electron ^-factor [51]. A related kinetic 
equation has also been applied in the context of thermonuclear fusion [[521 . In this model, however, both 
the dipole and the quantum dipole terms are missing. This equation can formally be retained from our 
result above by neglecting all terms of order h or higher. 



3.3 Thermal equilibrium 

The scalar spin distribution does also have some other important differences compared to the distribution 
for a classical dipole. An important example is the thermal equilibrium distribution. For in the classical 
case it would be given by 



/d(v,s) =nof M {v) 



1 HbBq 
4tz k B T 



sinh 



/ fi B B \ 
\k B T ) 



i -i 



exp 



k B T 



(58) 



where iiq is the equilibrium density, /m is a normalized Maxwellian velocity distribution, k B is Boltzmann's 
constant, T is the temperature and Bo is the external magnetic field. This distribution gives rise to a 
magnetization 

.T „ f HbBq^ 



V k B T 



(59) 



where 77 is the Langevin function. In the quantum mechanical case the corresponding distribution function 
becomes (assuming that the chemical potential is sufficiently large so that Landau quantization can be 
neglected, see Ref. |49l ) 



1 

Alt 



1+tanh (^V os0s 

\k B T J 



(60) 



where /m it the usual Maxwell-Boltzmann distribution. The zeroth order magnetization in this case is given 
by 



M, 



3ji B / dQ.sf q = «ojUfitanh 



V k B T J ' 



(61) 



as expected. The dynamics of the magnetization can be treated by the use of a fluid moment hierarchy and 
we now go on to consider this. 



4 Spin Fluid Moments 

Many problems do not require the full machinery of the kinetic approach and calculations can be greatly 
simplified if executed in a macroscopic fluid model instead (see for example [|53l and references therein). 
In particular when dealing with nonlinear problems the kinetic approach soon becomes very cumbersome 
[refs?] and the need for a simplified theory cannot be understated. 

In this section such a theory will be presented, derived from the kinetic theory by taking moments of 
the quantum kinetic equation l49l in a way analogous to what is done in classical plasma theory (see for 
example l53l ). The theory presented here was derived in l57l . As in the classical approach all intrinsically 
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kinetic features such as Landau damping will be lost, and one is faced with a closure problem, since the 
fluid hierarchy is an infinite series of equations that needs to be truncated at some point. 

Quantum fluid moments derived from the Wigner formalism have been applied before, see for example 
Refs. I54ll55l . and also in the case of gauge invariant Wigner-Stratonovich formalism l56l . In the spin-1/2 
case the moments have also been calculated ll59l starting from a quantum kinetic equation [[58] . In their 
treatment they use a matrix form of the kinetic equation and they also retain the collision terms which we 
here will neglect. However, their fluid hierarchy is only discussed shortly and also it is derived to order h 
and hence, for example the effect of the dipole term on the dynamics was not retained. 

Since we are working with a quasi distribution function in a phase space extended to also include 
the microscopic spin variable, the classical approach must be slightly modified. Firstly we also need to 
integrate over the microscopic spin variable, and furthermore a new macroscopic spin variable is defined, 
leading to a new hierarchy of spin-dependent macroscopic objects l57l . Thus we define the moments as 



dSlf, (62) 

- / dOxf, (63) 
n J 

m J dQ.{vj- iii) (vj - u j)f, (64) 

Qijk = mjdQ.(vj- Uj) (vj -uj) (v k - u k )f (65) 

S = - fdCl&f, (66) 



Pij 



n , 

Lij = ml da(3si-S t ){vj-Uj)f, (67) 



Ay* = m J dQ.(si-Si) (vj -uj)(v k - u k )f. (68) 

Here the first four moments are respectively the density, the fluid velocity, the pressure density and the 
energy flux density. The equation moment above d66l ) defines the spin density S = S(x,f) which yields the 
average spin density at position x and time t . The factor 3 in this definition occurs due to the correspondence 
(148V The sixth moment (l67l i is a mixed moment of the velocity and the spin which will act as some kind of 
spin-pressure. Finally, we have a mixed spin-velocity-velocity moment which could perhaps be termed the 
spin-pressure correlation. Similarly we could go on to define even higher order moments like the energy 
flux or a higher order mixed moment. Note that there is no need to include higher order moments in the 
spin variable since we have that J dClSiSj °< 5y. 

Using the evolution equation for the extended phase space distribution function $55[ we may now 
calculate the evolution equation of the different moments 



-nV-u, (69) 



dn 
dt 

du: q , U dB; 1 dPj; 

— - = ^(E i + £i j kU jB k ) + ^ Sj ^ (70) 

dt m m axi nm axj 

dPj; dU; dUi du k q q u dB k 
1 — —P., J - — Pr - — P- i + _p- P R +—F- P R 4- — £ - 

dt ax k ax k ax k m m m oxj 

M y dB k dQijk 

~/'*"5 5 ' ( /i - ) 

OXj ox k 



in 



dSi 2jU 1 dXij 

-j- = -i~ £ ijkSjB k ^ — (72) 

dt n nm oxj 

dY "' E dUk E dUj P dSi + V F y R 2fl R+Un dB < 

—r- — ~ L il^\ L ik^\ Hk-z: 1 £j k iLj k Di + — £,A/2^;#/ + jln- — 

dt ax k ax k ax k m n oxj 

dB k dAjfk 

-H nSi s k -±-—V±, (73) 
oxj ax k 
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where d jdt = d t + v • V. It is possible to derive the evolution equation for the Qijk, see ll55l and Ay^, 
however, these moments will then couple to even higher order moments. Instead some approximation is 
needed. In our paper ll57l we used the closure Qjjk = and Am = 0, which is perhaps the easiest one to 
deal with. 

As can be seen in Eq. d72i i. the spin-velocity moment E i; - acts as a pressure term in the evolution 
equation for the spin density. By including this moment, it is possible to capture some kinetic effects in 
a fluid theory. Some problems which might be difficult or tedious to solve in a kinetic theory may be 
reachable within a fluid theory. However, the exact role of the spin-velocity moment is a subject of further 
research. 

4.1 Two-fluid model 

It is shown in the previous section treating the kinetic equation that the distribution function can be divided 
in two parts, one for each spin direction along the magnetic field. Each fluid is seen to obey the same 
hierarchy as above f57l . and we will just have two sets of these equations, one for each species. Of course 
we will have separate macroscopic quantities n a , v a , S a and E a for each species. Here the subscript a 
indicates which species the quantities refer to. This approach adds a bit of complexity but can capture 
some kinetic effects due to a the different dynamics of the two spin states, and is therefore worth pursuing 
in problems were such physics is expected to play a role l68l . 

5 Applications of fluid and kinetic models 

The various quantum models developed cover several physical effects. Effects of the Fermi pressure and 
particle dispersion has been described in some detail in Refs. lfl"2"l [Tl both within fluid theories and within 
a kinetic approach. The fluid approach uses the so called Bohm-de Broglie potential to get an effective 
quantum force in the momentum equation, and the equation of state is chosen such as to get the Fermi 
pressure. Several modifications of classical behavior due to such models has been described in the litera- 
ture, see the references in Ref. Q for an up-to-date list. Many papers using kinetic approaches cover the 
effects of particle dispersion and Fermi pressure, but using a Wigner function derived without the magnetic 
dipole coupling of the Pauli Hamiltonian. In this way all effects due to the magnetic dipole force and spin 
magnetization is left out. In contrast to the these works, we will here focus on physical effects directly 
associated with the spin-coupling in the Pauli-Hamiltonian, which give raise to the magnetic dipole force, 
the spin precession and the spin magnetization in the above presented models. Most of the recent results 
along these lines has been derived from models similar to those presented here. 

5.1 Results from spin-fluid theories 

The most basic question to ask concerning the electron spin-properties in a plasma is "when are they im- 
portant?". For spin effects due to the Fermi-pressure this is straightforward to answer. The Fermi -pressure 
becomes important when the Fermi temperature approaches the thermodynamic temperature, which give a 
simple condition on the temperature and density of a plasma. For the effects due to the direct spin-coupling 
in the Pauli-Hamiltonian the answer is less straightforward, as it depends on the full parameter regime 
(involving also the magnetic field strength) but also on the specific geometry of the fields. During certain 
geometric configurations, the spin effects can be important in regimes of modest density and modest tem- 
perature, which traditionally has been thought to be completely classical l60l . A specific example of this 
kind can be found in the MHD-regime. In Ref. [60] fluid equations of the type d25ll-(F26]l was adopted to 
MHD-regime, in order to study the physics of nonlinear spin-modified Alfven waves. Within linear theory, 
the Alfven waves was almost unaffected by the spin terms, provided the Zeeman energy associated with the 
unperturbed magnetic field Bq was much smaller than the thermal energy, i.e. /ig#o *C kgT . This condition 
holds for most plasmas except close to pulsars and/or magnetars. However, nonlinearly the situation is 
different. Ref. [60] used a two-fluid spin-model based on d25ll-(l2"6ll, where spin-up and spin-down popula- 
tions were formally treated as different species, as a means to capture certain kinetic effects within a more 
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Figure 1: A figure illustrating regions of importance in parameter space for various quantum plasma ef- 
fects. The lines are defined by different dimensionless quantum parameters being equal to unity. The effects 
included in the figure are: (i) Fermi pressure effects, described by the parameter Tf /T e °^ h 2 n^/ 3 /inksT,,; 
(ii) Effects due to the Bohm-de Broglie potential, described by the parameter hco pe /kBT e ; (Hi) Spin single- 
fluid Alfvenic effects, described by the parameter h 2 (0 2 /mc 2 kBT e ; (iv) Spin single fluid acoustic effects, 
described by the parameter [IbBq / 'ksTf.. Here lines for three different magnetic field strengths are drawn. 
The quantum regime correspond to lower temperatures, i.e. it exists below each of the three horizontal 
lines; (v) Spin two-fluid nonlinear effects, described by the parameter HsBo/miC^. The quantum regime 
corresponds to higher densities, i.e. it exists to the right of each of the three vertical lines. 



simple fluid theory. From this theory a nonlinear Schrtidinger equation 

v' , IB, I 2 

id,Bi + -±d 2 B x + Q^-Bi = (74) 
I B Q 

for Alfven waves propagating parallel to the magnetic field was derived, where B\ is the slowly varying 
magnetic field amplitude, v' = dv g /dk is the group dispersion, £ = z — v g t is the comoving coordinate and 
v g is the group velocity. These quantities are determined from the Alfven wave dispersion relation, which 
reads a> 2 = k 2 c\{\ ±kcA/(O c i), when weakly dispersive effects due to the Hall current is included |6T1 . 
Here is the Alfven velocity and co c , is the ion-cyclotron frequency. The upper (lower) sign corresponds 
to right (left) hand circular polarization. The nonlinear coefficient is Q = Q c [\ — / ' m i c l) 2 ]> where 
the classical coefficient is Q c = kc\/A(c\ — c 2 ) ~ — kc\/Ac 2 , where c s is the ion-sound speed. Although 
linearly the modification of the Alfven waves can be neglected (for modest temperature and densities), 
the nonlinear coefficient Q could be significantly affected by the spin terms. Illustration of the parameters 
needed to make the different quantum plasma effects significant, is shown in Fig. 1. In particular we note 
that the two-fluid nonlinear spin effects are important for high plasma densities and/or a weak (external) 
magnetic fields. For comparison, both the Fermi pressure and the Bohm-de Broglie potential need a low- 
temperature or a very high density to be significant.. A somewhat surprising result is that here nonlinear 
spin effects tend to be more important for a lower magnetic field, whereas the opposite is true for linear 
spin effects. 

As a second example of spin fluid effects we will consider the ponderomotive force of an electromag- 
netic wave propagating along an external magnetic field . Classically, the density fluctuations induced by 
the ponderomotive force of an electromagnetic (EM) wave lead to an electrostatic wake field l62l . as used 
in advanced particle accelerator schemes |63l . In other regimes, the back-reaction on the EM-wave due 
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to the density fluctuations leads to phenomena such as soliton formation, self-focusing or wave collapse 
[64 65]. When spin effects based on Eqs. (|25]l-(|26]l are included, the spin-contribution to the pondero- 
motive force, resulting from the combined effect of the magnetic dipole force and spin-precession, leads 
to a separation of spin-up and -down populations. The expression for the ponderomotive force density for 
electromagnetic waves propagating along the unperturbed magnetic field Bo = Bqz can be divided into its 
classical part 



Fez 



2m 2 (0 ((0 ± (0 C ) 



■± 



k(Q c 



dz 0){(0±(Oc)dt 



where +(— ) correspond to right (left) hand circular polarization, and the spin part 



Faz = T 



2fl 2 Spg 
mh (co±co g ) 



d_ 



((O±C0g) dt 



\B\ 



(75) 



(76) 



Here E and B are the electric and magnetic field amplitude, respectively, at c = eB^/m is the cyclotron 
frequency, co g = 2jiBo/h = g(8 c /2. The index a = (u,d) refers to the up- and down populations, respec- 
tively. In particular the up- and down spin vector components are So,, = 1, and Sq c i = — 1. Thus it is clear 
that the spin-part of the ponderomotive force induces a separation of the up- and down-populations. This 
effect survives also in the absence of an external magnetic field, and it turns out that the magnitude of 
the relative up-and down density perturbations can be larger than the classical density perturbations in an 
unmagnetized plasma, provided 



1 < 



hco (0 2 p L 2 



(77) 



mc^ c 

where L is the pulse length of the high-frequency pulse. The first factor of the right hand side of (ITTl i is 
smaller than unity for frequencies below the Compton frequency, but the second factor can be large for 
long pulse lengths, and hence large spin-polarization can be induced by a sufficiently long EM-pulse in an 
unmagnetized plasma for optical frequencies and higher. Once the plasma is spin-polarized, the spin terms 
in the evolution equations can be important for the dynamics. For further studies of results from spin-fluid 
theories, see e.g. Refs. 



5.2 Results from spin-kinetic theories 

The full kinetic theory dBTl i is accurate, but cumbersome for many purposes. As seen from the full (15 It . 
the effects separates quite naturally into particle dispersive effects (which are insignificant for spatial scale 
lengths much longer than the characteristic de Broglie wavelength), and effects due to the electron spin. The 
particle dispersive effects has been studied in some detail in e.g. Refs. lfT2l [T1l. Focusing on the spin effects 
rather than particle effects, we can therefore consider the long scale-length equation (|55l . An interesting 
effect of spin-kinetic theory is the appearance of new wave-particle resonances. Linearized theory in a 
magnetized plasma can be solved in much the same way as in a classical plasma (see e.g. Ref. lIBTl for 
technical details). However, due to the fact that the spin-precession frequency co g and the Larmor gyration 
occurs (£> c are slightly different (i.e. (0 g = (g/2)co c « 1.001o c ) new resonances appear. In particular the 
denominators of of the perturbed distribution function in kinetic theory are replaced as 

7 — v ^ t — ; v < 78 > 

(o - k z v z - nCO c ) ((0 - k z v z - nca c - mco g ) 

when spin-kinetic effects are included, where the integer n covers ±°° and m = ±1. Thus for a fixed 
parallel phase velocity, resonant wave-particle interaction can occur at a much lower temperature when 
spin effects are taken into account, since \(O c — C0 g \ |o) c |. This aspect has been discussed in some detail 
in Ref. J49). Furthermore, for perpendicular propagation to the external magnetic field, new Bernstein- 
like modes appear with frequencies close to the resonant value CO = Ao) ( = | (O c — (0 g | . Specifically the 
dispersion relation for perpendicular propagation, reads 



co 2 = k 2 c 2 + co 



I j {j 2 (k ±V J(O c ) 



k z n z A(O c sin z 9 s 
4m e (co - Aco c )k B T 



[J 2 (k ± vj(0 c )]}foda. (79) 
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where Jo and J\ are zero and first order Bessel functions, respectively. Here we have assumed the clas- 
sical terms involving higher order Bessel-functions are negligible, which is accurate for (0 <C | C0 C \. The 
numerical solutions of (|79l reveals that the wave frequency only deviates slightly from the resonance Aco c 
for most parameters J3T). It should be noted that a Madelung approach (c.f. Eqs (f25ll-(t26l>^> cannot capture 
the physics of the resonances at (0 w Ao c . However, the moment theory (|69l)-(l73b correctly recovers the 
dispersion relation in the low temperature limit 1571 . Due to the higher order moments, however, it should 
be noted that this theory is computationally somewhat more demanding than Eqs. (|251)-(|26T). at least if the 
higher order spin effects are omitted in that theory. Other works on spin kinetic theory includes Ref. |70l . 
where the general linear theory based on (l55l > was studied in a magnetized plasma, Refs. Il52l . where a 
simpler version of (1531 ) (a semiclassical correspondence without the magnetic dipole force term) was stud- 
ied with regard to fusion applications, and Ref. [69] where a semiclassical version of d55l l was adopted to 
consider spin induced damping of electron plasma oscillations. 

6 Summary and Discussion 

A rapid development of spin-models for plasmas has taken place during the last few years, using different 
types of methods. The most accurate are the kinetic approach, based on the combined Wigner and q- 
transforms of the density matrix. This theory captures not only the spin effects, but also the particle 
dispersive effects. However, for scale lengths longer than the characteristic de Broglie wavelengths only 
the spin effects remain, as described by d55T >. While the kinetic approach is accurate, and contains much 
interesting physics, as for example the resonances displayed in Eq. (l78l l. there is a simultaneous need for 
models that are simple enough to be applied also for more complicated problems, involving inhomogenities 
and nonlinearities. Various fluid approaches have been adopted for this purpose. Those based on the 
Madelung approach capture most of the spin-physics in the physically intuitive effects of a magnetic dipole 
force and spin precession, together with a spin magnetization current. Such approaches has been used 
in several recent works, see e.g. Il66l l67l l68l . However, higher order quantum terms also appear in this 
context, which either must be modelled or omitted, in order to form a closed set for the fluid variables. 
Another approach to obtain fluid theories is by computing moments of the kinetic theory. The basic terms 
(i.e. spin-precession and the magnetic dipole force) are the same as in the Madelung approach, but now a 
spin-velocity correlation tensor appears in the evolution equation for the spin. This has a correspondence in 
the tensor of (l28l l in the Madelung approach, but here things are complicated by the presence of other terms 
such as given by d29l . Thus an advantage with the moment approach is that it is straightforward to model 
the spin-velocity tensor by computing the corresponding moment. As always when taking moments, the 
coupling to a higher moment appear, but truncating the moment expansion in the next step (i.e. dropping 
higher order moments in the evolution equation for the spin-velocity tensor) seem to capture most basic 
spin physics, and leave out only thermal effects l57l . However, even when one is interested in the low 
temperature limit, certain effects of a nonzero temperature should be kept to address the behavior in the 
more common plasma regimes. This has to do with the fact that the Zeeman energy is typically much 
smaller than the thermal energy (i.e. HbBq <§C kgT), such that in thermodynamic equilibrium the two spin 
states are almost equally populated. As a consequence, even if the temperature is small in all other respects 
(i.e. all characteristic velocities of a system is larger than the thermal (or Fermi) velocity), there is a 
large spread in the spin distribution. As a means to capture some of the physics associated with this large 
spread, two-fluid theories of electrons has been developed where up- and down states with respect to the 
(unperturbed) magnetic field are described as different species, as described briefly in section 3.1. To some 
extent the fluid theories including the spin-velocity tensor seem able to account for some of the effects of 
the spread in the spin distribution, but it is too early to make a definitive evaluation of the various models 
strengths and weaknesses. Thus the final conclusion is that more research on the physics of electron spin 
in plasmas is needed. 
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